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15. We shall now apply the equations of Article (12) to 
the Lunar Theory, in the following manner. We shall first 
suppose that the moon moves in the plane of the ecliptic, and 
afterwards prove the results obtained on this hypothesis to be 
true when the inclination 1 is taken into account; we shall then 
determine the motion of the plane of the orbit, the mean 
motions of the perigee and node to the third order, and the 
parallactic equality. | 


By Article (12), we have a 
fr 


| dh dR 
and by Article (3), we have 
dt? dé dt 


‘Therefore (1) becomes 
| d*u u r aR 1drdR 


and since = (2) becomes 
d (h’) dk 4 
16. To determine R, we have, by Article (14), pons 
cos (7, 7”) = p, 1 
R=m 
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2 


(1 — 2pz +2") where z = 
r 


l 


_ Hence, retaining only the lowest power of z or a9 We have 


(observing that p or cos (7, 7’) = cos (0 - 8’), since we suppose 
the moon to move in the plane of the ecliptic) ) 
m 1 +3 cos 2 (0- @) 


4 
my 3 sin 2(0- 60) (6), 


We shall expand ‘the second members of (3) and.(4) as far 
as the third order of small quantities, retaining, however, only 
those terms of the third order in which the coefficient of @ is 


nearly equal to unity or zero ; and this we shall do by substi- © 


tuting for 7, 7’, and @' their values in terms of @ to the first 
order, namely | 


=m9+ 2e sin (m0 + ...... (8). 
17. To determine 7 in this manner, we have, by (4) and (6), 
ain 3 (0.- @); 
or, y ( pu 5 | 
{1 — 4e cos (8 —w) + cos (6 —w')} sin 2(0- 6’). 


In the terms of this, which are multiplied by e and e’, we 
may put for @' its first approximate value, namely, mO + p. 
In the remaining term, if we put for 6’ its value (8), expand 
to the first power of e, and then transform products of sines | 
and cosines into sums, in the usual way, it is evident that 
sin 2 (@ — 0’) will contain no term in which the coefficient of 
§ is nearly equal to unity or zero. It follows therefore that 
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“we may suppose 6’ to be equal to m@ + B, simply, i in all the © 
terms of this equation. 

This being the case, the product cos (6 — w’) sin 2 (0 - 9’), 
transformed into a sum, gives no term to be retained, and the 
product cos (@ w) sin 2 (0- @) pres only } sin (6 - + w): 

_ we have therefore 


= {sin 2 (0. 9) - 2e sin 26 + 


Therefore, integrating, ‘supplying the proper constant 
{au (1 e*)}, and neglecting powers and products of m and 
above the second order, we have 


= ap {1 4 2 (0 - 6m’e cos (0-20 


= cos 2(0- 6m*e cos (0 201 + 


18. To find - . , we have, (5) and 


| 
= {1-Secos(O- w)+ w')}. ) 


7? Zor exactly ie same reasons as before, we may here suppose 
# to be equal to m@ + 3 — and then, retaining only the 


terms, we have 
= {- cos 2 (0 - 6') +e cos(0—w) 
= + 2ecos (0 — 20 +) cos (0 - w)}. 


19. Lastly, to find 7, tS ,we have, by (6) and (7), pro- 


' ceeding exactly as before, 
1 dr dR. 3m'r’ dr | 
h? dd apr” dd ‘ 
—w) sin 2 (6 - 0’) 
20 


€ COs (0 - 
4a | 


20. Let us now, for brevity, put a 
20-0)=9, 0-20+=y¥, (9). 
here being simply + 


a oz 
a 

4 
: 
| 
| 
| 
» 
2 
§ | 
| 
| 
wet 
| 

: 
et 

as 
4 
2s, 


~~ On the Lunar Theory. 


Also, for e cos (8 - w) put aw- 1; and then we have 


1+6. 


3 


om { 2-5 cos cosy 5 cos x| 
h* dd 2 

1 dR wv { 3: 


Hence, (3) becomes 


(, _ 
| 


= ; (1 + 2m? 3m’ cos + mie cos 3 me’ cos X:) 
Assume now putting 1 - = 


fe cos - w) + A+ B cos Coos y+ D cos} 


hee, by (9), 
ou. 
+ _ 2m)} B cos 


Hence, equating coefficients, neglecting powers and pro- - 


ducts of m and e above the second order, and putting 


(1+ 2m’) =u,, we have 


U=U, +e cos (cB - = cos 
21. To find ¢ in terms of 0, we have 
| dr 


in which we must substitute the values of A and r already 
found in Arts. (17) and (20), and then integrate. But since 
h is found by an integration with respect to 6, and ¢ by 
another integration with respect to 0, it is not sufficient to 


retain merely the terms of the third order in which the 


coefficient of @ is small in the expression for — # , but we 
must also preserve all similar terms of the Jourth order.* 


* In Airy’s Tracts and Pratt’s Mechanical Philosophy the terms of the third 


order only are retained, and therefore it appears to me that ¢ is found by an 


inconclusive method, so far as terms of the second order are concerned. 
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To obtain these terms in the expression for a , we must 
substitute for 7*, 7° and 0 their values to the second order, 


namely 


{1-8e cos (6 


— 
_ (we may look upon e” as of the same order as m’). 
+ B+ sin (m0 + B ~ - 2me sin (8 w). 


Now it is easy to see that, on account of the expression for 


ra being multiplied by sin 2(6 - @), the expansions of 


rand @' will give rise to no terms in which the coefficient of 


6 is small; we may therefore suppose r =a and 0 =m0+ 8B 
simply. But in the expansion of r*, the term 5e* cos 2(6 - w) 
combined with sin 2 (0 - 6’) gives a term to be retained, 
namely 3 e sin We have therefore 


+B me’ cos 2(0' - w)}. 


Introducing this additional term in Art. (20), we find 
u=u,{...... cos 2 w)}. 

22. We may now proceed to find ¢. 

Represent the values of w and kh’ by u, (1+ U) and 
- e)(1+ H); then we have, by (10), 
1 
ao” u,v 
The values of U and H (neglecting in H the term of the 
third order in which the coefficient of @ is not small) are 


U = e cos — w) + m’ cos 


20! 
cos C08 x + C08 2(0' - w), 


cos — cos 2(0'- 6); 
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U* evidently contains one term of the second order and two 


of the third, 
| cos’ (c8-w), 2m'e cos w) cos ¢, 


and mee cos (cO — w) cos y. 


The second of these gives no term to be retained, the third 
gives (putting = 1 and observing that = - 20° + w), 
— cos 2 (6 - " and the first is equal to 


{1 + cos w)}. 


Hence the of the term cos 2(6' - 0 


(11) is 3+ 4, which = 0; this term therefore disappears. 
Therefore, collecting the terms of (11), and putting all the 
] 
constant part = —, we have 
n 
ndt 


1 - (cB - w) 
4 11m’ 


oe lome 
+ 008 2 (8 w) - 


cos + 8m’e cos x; 
and therefore, integrating, ; 
nt + const. = 2e sin — w) 


+ sin ~ sin - sin 8me sIn w. 


23. We now proceed to shew that these results are not 
affected by the inclination of the moon’s orbit to the plane of 
the ecliptic. 


By. Art. (14) p is the only part of # which contains :; 
therefore, since, by Art. (16), | 


2 
we have (12) 


| and, by Article (14), observing that \’ = 0, 0, = 0, 
p = cos cos (0, - &), 


sin cos (6, - — - cos sin (6, .(18). 


Now, by Article (14), - 


tan (6, - ») = cos. tan (0 »), 


4 
a4 
% 
Va 
& 
ae = 
A 
Ross 
4 
a 
ve 5 
4 
ae 
« 
a 
2 
“ed 
| 


On the Lunar Theory. 


hence, differentiating and then neglecting higher powers of 


and 1, and putting 6, = 0, we find 


sin sin’ (6 — cos =- sin (6 - v) cos (0- 
t t 


therefore, by (13) and (12), 


| v) {sin (0-v) cos (0-@) - con sin 


dR 38n'r’ 


a r 


Hence, if we heh by R, the part of R independent of the 
inclination, we pitt integrating, 


* cos (0 sin (0 -») 


cos (0 sin (0 - y) sin (0 - (14). 


from which it appears that the part of R depending on the 
inclination is of the fourth order ; therefore the only part of 
the preceding investigation that can be affected by the in- 


clination is the determination of I? from ae ‘Article ( 17), 


- and we need only look out for. those terms in which the 


coefficient of # is small in = Now 
cos (20-6 - v) sin (0 


The part multiplied by @ here evidently gives no term in 


which the coefficient of @ is small, we have therefore _ 


dd 
and therefore the value of h’, found in Article (17), is the 
correct value when the inclination is taken into account. 


Again, in finding the time, instead of h=?r° 4 , we have, 


by Article (12), 


= 


le (12 
by Anticle ). 
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The part involving « here is of the fourth order, and is 
therefore to be neglected altogether. 

Hence it appears that the supposition that the inclination is - 
zero does not in any way affect the process of finding the 
values of « and ¢ in terms of @ to the second order of small 
quantities. It is important to remark that @ is the longitude 
of the moon measured on the plane of her orbit ; and to in- 
troduce the longitude measured on the plane of the ecliptic, 
namely 0,, we must apply _ common reduction 


5 sin 2(0-»). 


This explains the reason why the only term in the expression 
for the time in Airy’s Tracts, which depends on the inclina- 
tion, is the reduction just put down. 


24. We now proceed to determine the motion of the 
plane of the moon’s orbit. 


By Article (12)* we have | 
@r ., h. 


therefore, by Article (23), equation (12), we have 
= cos (0 - sin (0 - ») sin (6' - »), 


=— sin (9 - @’) sin (6 ») sin (0 - 


In both these put 
2 sin (0 — v) sin (0 v) = cos (6 - 8) cos (0 + 0 — 2y), 


retain only the terms which rise in integration, and the 


constant term in —, and we find 


dé 
{1 — cos’ (0 v)}, 
sin 2 (0 v), 


3m 
which give = const. 0+ + sin 2 (' »), 


= const. + 008 2(0' - 


* Ih Arts. (12), (8) and (9), in the expressions for © —, Sintis put by mistake | 
for tan t. 
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and these values, substituted in the equation 
sin A= sin, sin (0 - »), 
give the moon’s latitude. | | 
25. We shall conclude. by pointing out briefly how the 


-parallactic inequality and the more correct values of the mean ~ 
motions of the perigee and node are obtained. © 


In Art. (16) the term of R, involving the third power of.2 


r 
or is 5p” - 


In this term neglect the eccentricities altogether, put for 
its value 3 cos 3 (0 0) + 3 cos (0 - @), retain only 


the part multiplied by cos (0 - 6), which rises in integration, 
and, so far as this term of R is concerned, we find | 


{8 sin 0-09}: 
hence in 73 there arises the term 
=, cos (0 (B); 


dr | 
and, observing that e is zero and therefore @ = 0, in the 


equation (3) there arises the term | dO 
(A) +(B) or (0 6), 


If, therefore, F' cos (0 — 6’) be the corresponding term in 4, 
we have 


and therefore | '= 
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Hence we have 
u= m cos (0 - 
16 


The term here obtained is the parallactic inequality. (See 
Airy’s Tracts.) 


26. To obtain the more correct value of the mean motion 
of the perigee, we must look out in the second member of 
equation (3) for terms of the form Ade cos (0 - w). Now 
it is easy to see that, when we substitute for 7, 7’ and # 
their values to the second order (see Art. 21), in order to 
expand the second member of (3) to the fourth order, » 
terms of the form Ae cos (9 - w) will arise only from the 
combination of the quantities me cos (or sin) (8 — 20 + w) and 
cos (or sin) 2(0 ~ @). Hence, confining our attention only to 
the term me cos (0 ~ 26’ + w), ¢.¢. putting for 7 


a {1 (0-20 + 


and for and a and wed we find (see. 
_Art..17)- 


hence, and putting au -— 1 for e cos (0 — w), we find 
we have (see Art. 18) | 
rdR m 15me 3 cos 2(0-@) 
ap 3 cos (9 20 + | 
m* 9.15 
9.15 
Lastly, we have (see Art. 19) 
dr 
= @ in (0 - 26 + 
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1 dr dR 3m 15me 
and (0 20° + w)sin 2(0- @) 
$3.15 


225 


4 8 8 16 
Hence the differential equation for «(see Art. 20) becomes 
du 3m? 225 ; 
3m" 225 
Hence we have ——-—— 
2 2,16 
3 2 
4 32 | 


| ore determines the inore accurate value of the niean motion 
of the perigee. (See Airy’s Tracts.) 


27. The more accurate value of the mean motion of the 
node i is easily obtained from the equation (see Art. 24) 


} 


by putting for v, in the second member, its first approximate 
value 


| om. 
+ sin 2 (0 - »,), 


which, expanding cos 2 — ), and only the con- 
stant part, gives 


Om 


oa determines the more accurate value of the mean motion 
of the node. (See Airy’s Tracts.) M. O. B. 


Hence, adding (4), (B) and (C), we have a term : 
— 
x 


II.—ON THE EQUATION (D+ a)" y= X. 


Tuts short paper is intended to illustrate the formation of all 
the differential equations connected with an equation of the 
nm order, namely, its successive primitives of the 
(n — 2), &c. orders. One of the greatest encouragements in 
the use of the calculus of operations is the perfect view which 
it gives of all the primitives of a linear equation of the form 
(D-a)(D-5b).... y = X, allowing any one of them to be 
singled out and formed independently of the rest. This is 
not so easily done when two or more of the set a, 6, &c. are 
equal, the reason being that some of the primitives in such a 
case are not linear with constant coefficients. 1 do not intend 
to consider all the combinations of primitives which may occur 
in an equation of the general form (D+a)"(D+0b/...... 
- y = X, but only the particular form (D+ a)" y = X. 

To take first the simplest form, let a = 0, or y) = X. The 
discussion of this easy case at length will allow of brevity in. 
treating the more general case. First, it is desirable to express 
in terms of {x"pzdz: let last be we have then 
the following theorems : 


r + 1) fdry" = ne P. + 


= B- B+ A'fdzy B- &c. 
From the last theorem (which is an obvious extension of 
John Bernoulli’s) it follows that if A be a rational and integral 
function of the n“ degree, {ABdz is always integrable when 
B can be integrated n+ 1 times. Hence pz .y is integrable 
per se, if pr its integral being 
px.y” &c. This must be called a known theorem ;* 
but it certaimly is not well known. 

The mode of proceeding with y” =X is now obvious 
enough: let our instance be and let P.. 
Multiply successively by 2’, x’, x’, x’, and integrate ; we thus — 
get the four primitives of the third ‘degree, which are (the 
constants being in the 4 Py &¢.) 


(1) Ps 

(2) ty” y = 

(3) — + =P,, 
(4) xy + - by = P.. 


* It is a particular case, deducible from the general method. of finding 
integrating factors for linear expressions. 
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To find the six primitives of the second order, observe that 
the first side of (1) 3 is integrable, and remains so when multi- 
plied by 2 or 2’; also that the first side of (2), or the same 
multiplied by z, is integrable ; also that the first side of (3) is 


integrable. Here are then the processes for finding the six 
equations required, which are 


(5) y =2P,-P,, 

(7) vy — 2ry' + 2y = P,), 

(9) Bey! + =] - P), 
(10) 2y' - 4zy' + 6y = 2P, - P,. 

‘To find the four primitives of the first order, observe - 
that (5) is integrable, and remains so when multiplied by 2, 
and that (6) and (8) are integrable. ‘The four equations _ 
required are 

(11) = 3(2'P, -22P,+P,), . 

(12) zy -y=52°P,- +! 

(13) 2y' 2y =) 2°P,-\ cP, +P, 

(14) - 8y=42°P,-2P,+} P,. | 

Finally, — of (11) gives the primitive of the order 


zero, 
(15) y= - 32°P. + 8aP,- 
To extend this to the case of oe + a)'y = X. Search in the 


common way for the factor which renders (D + a)'y integrable - 
per se; it will be found to be e“(¢,+¢7+.... +6, ,0""), 
say «“. pz. ‘To find D". «“px (D+ a)" y, transform it into 
Dy or D" (pxD"e“y), which, on account of the 
degree of pz, appears in finite form thus, 

Proceed thus with (D+ ay y= X. Let be 
by P.. We have then, as 


mintl 


Ja" = 


(D+a)". AB= A(D +a)'B-A' (D + 
Though it may be worth while to exhibit these theorems in 
terms of the operation D+ a, yet it 1s not necessary to pro- 
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further ; for (D + ay y= X is = X, and the 


rule of proceeding is to solve y) =e" X, as in the manner 
above given, and then to write in the several equations ey, 
+a)y; «”“(D+ayy, instead of y, y', y’', &c. Thus, 
~2y' = Xaxdz { Xx'dx being one of the of 
y” = X, it follows that 


(D+ ay y — (D+ a) y= [Xe“z'dz, 
or xy'+(2az- 2) y'+(@a- 2a) y=e wal Xe “edz — 
is one of the preeoei of 
+ + + 4a’y! + = X. 
The following theorem corresponds to that of integration by 
parts + a)’ {A (D+ a) B} = AB-(D+a)' (BDA), 


which is an immediate consequence of the simple 5 acai 
relation D 


A. D. M. 


III.—ELEMENTARY DEMONSTRATION OF DUPIN’S THEOREM. 


Ir there be three series of surfaces, such that all the surfaces 
of each series cut the surfaces of the other two series at right 
angles, the theorem to be proved 1s that the lines of inter- 
section of any one of the surfaces of the three series, with the 
surfaces of the two conjugate series, are its lines of curvature. 

Let O be any point in which three conjugate surfaces 
intersect, and let the rectangular axes OX, OY, OZ be 
: perpendicular to the tangent planes of the three surfaces 


at O. Let . (a) 


Fn (Bs Ys = (@) 


be the equations of the three series ; and, when proper values 
are attached to A, A,, A,, let (a) be the surface touched by 
YOZ, (a,) by ZOX, ‘and (a,) by XOY. 


Hence, when z=0, y=0, z=0, we have 


dy dz 

GY, , 
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= Now, since the system is orthogonal, — 
GF, UY, LY, _, 
"dy dy dz dz. 
dz dy dy" dz dz 
fd, 
dx dz "dy dy” dz dz 


‘Differentisting the first of these equations with respect to z, 
the second with respect to Y> and the third with respect to z ; 


putting z, y, and z each = 0, in the result, and making use of 
equations (b), we have 


the brackets denoting that in the quantities enclosed, z, y, z 
are put=0. From these equations we conclude that 


are each = 0. 


Now, since YOZ is the tangent plane to the surface (a) at 


=0 


O, we hake 
d*x d’x\ .\ 
for points in the surface adjacent to O. To determine (5 Z| 
we have, from (a), | dydz 
ie dy 
dy of 


Hence, differentiating with regard to 2, and putting z, y, 
_ and z each = 0, in the result, 
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Hence (d) becomes 


Hence the planes of zy and zz contain the principal sections. 
of (a) through O, and therefore the lines of intersection of | 
(a,) and (@,) with (a), touch the principal sections of (a) at O. 


Now O may be any point in one of the surfaces (a), and 


therefore each of these surfaces has its lines of curvature 
traced upon it by the surfaces of the series (a,), (@,). Similarly 
- it may be shewn that each surface, (a,) has its lines of curva- 
ture traced by (@,) and (a), and each surface (a,) by (a) and 
(4,), which is the theorem to be proved. | 
sy differentiating the first of equations (c) with regard to 
y and z respectively, and putting z, y, and z each = 0 in 
the results, and performing corresponding operations on the 
second and third, we obtain expressions for 


which lead directly to the expression for the curvatures of the 
principal sections of the three surfaces at O, given by Lamé, 
in his Memow on Curvilinear Co-ordinates. 2s. 


IV.—NOTE ON A DEFINITE MULTIPLE INTEGRAL. 
In the XVIII number of this Journal Mr. Boole pointed — 


out the incorrectness of a theorem given by M. Catalan. The 
following pages contain a brief demonstration of the result to 
which he was led. Both he and M. Catalan made use of 
what is generally known as Liouville’s theorem, and thus 
perhaps rendered their analysis less simple than it would 
otherwise have been. 

Let us transform the integral 


by the assumption 


a2 +....4,£,=au, . (a = 3a’), 
u,...+.%,, and such that 
= Se. 
Then, as is well known, dz,.... dz, is to be rhea 9 


du,.... du,, and thus 
fdz,.. [du, fau, fdu,... fdu,...(1 
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Let the integrations on the first side of this equation include 
all values of thé variables which do not transgress the limits 


= 
where B is supposed to be greater than A. Then, as S2*= Su? ; 
the corresponding limits on the second side of the equation 


are | = = B’. 
Transform the integral in z by assuming 


z=rcos@,, z,=rsin 0, cos sin 


and that in w, by similar assumptions, —— 
u,=rcos9’, u,=rsin cos %,...u,=7 sin 6,...sin 


I may be allowed to mention that this tinnsteialiim, 
which appears to have been given for the first time by 
Mr. Boole, in the last number of the Journal, had occurred 
to me before I had seen his paper. His analysis leads at 
once to the conclusion that dz,...dz, is to be replaced by | 


‘This may be also proved by successive lide in the 
manner pointed out in the case of three variables by Mr. A. 
Smith, in the first volume of the Journal. 


Thus (1) becomes, since = Xu’ 
fsin””6, d6... a0, (a, cos 0, a, sin 6... sin oF 


_fsin 
The limits for @ and 6’ are ne same. 
That this equation may subsist for all values of A and B, 


it is necessary and sufficient that — 
... [d0,_, f{r(a, cos 0,+...a, sin 8, ... sin ,)} 
= f(ar cos 6',) d6’, fsin”*6',d0’, ... 10 bs 


With respect to the limits of @ and @’, it is i difficult to 
perceive that if 6... are taken the limits 0 and 
Z,...%,, will receive all the values of which they are capable, 
namely, all that included between - r and + 7; and that the 
same set of values cannot occur more than once. But in 
order that z, may vary from 0 to — 7, it is necessary to extend 
the superior limit of 0, from m to 2x. Thus the limits of 
are 0 and 27, while those of the other variables 9,...0, , 
are 0 and And similarly for @. 
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66. Note on a Definite Multiple Integral. 
~ On the second side of (2) we have the factor 


Let be odd, then we shall have o 


Sunilarly 
sin” db’ sin” 9 = &c. = &e. 
| 
Lastly & 
7 2 
Thus (4) = 2 


Again, if m be even, we get 


(A) = — 


Ana this, multipled by or gives, as before, 


(A) - 
and thus (2) becomes, making r=1, | 


sin””*6,d6, dé _, f(a, cos 0, + ... &e.) 
| | 


(AS) : 


2 


This is the general result given at the end of Mr. Boole’s 


paper, and which includes the others. 6 oe 
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1V.—NOTE ON SOME POINTS IN THE THEORY OF HEAT. 


In problems relative to the motion of heat in solid bodies, the 

initial distribution, which is entirely arbitrary, is usually one 
of the data. When this is the case, and the circumstances in 
which the body is placed are known, the distribution at any 
subsequent period is fully determined, and if our analysis had 
sufficient power, would become known in every case. ‘The 
solution of the problem would be an expression for the tem- 
perature of any point in the body in terms of the co-ordinates 
of the point, and the times measured from the instant at 
which the distribution is given. — 

It is in many cases an interesting investigation to examine | 
what this expression becomes when negative values are 
assigned to the time. If, for a particular negative value 
— r, we find that the expression gives an actual arithmetical 
value for the temperature of every point in the body, the 
distribution represented will be such that, if it had existed 
at a time +r before the initial instant, the given initial dis-— 
tribution would have been produced by the spontaneous 
motion of the heat. 

It is clear, however, that the arbitrary initial distribution 
may be of such a nature that it cannot be the natural result 
of any previous possible distribution, or that it cannot be any _ 
stage except the first in a system of varying temperatures. 
This, for instance, will be the case if there be any abrupt 

transitions in the initial temperatures of adjacent points, or if | 
the curves representing the initial temperatures of points 
‘situated along any straight line through the body, have cusps 
or angular points ; for, though we may suppose such a dis- 
tribution to be arbitrarily made, it could obviously not be 
produced by the spontaneous motion of heat from any other 
preceding distribution, and all such abrupt transitions or 
angles which may exist in an initial distribution, will dis- 
appear instantaneously after the motion has commenced. 

If, however, in any case, a complete solution of the problem 
has been obtained, that is, if an expression for the tempera-_ 
ture of any point, in terms of its co-ordinates, and the time 
has been found, this expression must assume some form, when 
negative values are given to the time; and it is my object 
here to examine the nature of this form, when the initial 
state is such as not to be deducible from a previous distribu- 
tion. | 

The simplest case is that of the linear motion of heat, and 
as from it we are enabled to understand the nature of the 
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problem in its most general form, I shall, for the present, 
confine myself to this case. The case of a thin rod, protected 
from any lateral radiation, or of an infinite solid, in ‘which the 
temperature is distributed in parallel isothermal planes, is 
what is usually contemplated in speaking of the linear motion | 
of heat; but the case of a thin rod losing heat by radiation 
from its sides may also be readily reduced to the same as the 
two former. For simplicity, however, we may suppose that 
there is no lateral radiation. 

Let a fixed point O in the rod be considered as origin, 
let the distance to any. point P in the rod be z, and let the 
temperature of P, at the time ¢ from the initial instant, be v. 
The equation of the motion of the heat will be 

do | | 
| 
if, ie brevity, we choose the unit of heat such that the con- 
ducting power of the body, referred to a unit of its volume, 
may be unity. 
The complete integral of this alia is 


Let ,v be the initial temperature of P. Then, putting 
¢= 0, in (2), we have. : 
(0 = ZA; cos + n,;) 
Now, whatever be the arbitrary nature of the function ,r, 
whether continuous or discontinuous, it has been shewn by 
Fourier that it may be represented by a series such as the 
~ second member of (8), and he has shewn how, when the 
value of wv, corresponding to any value of z, 1s given, the 
constants A;, m;, m; may be determined. ‘The series thus 
found is an actual quantitative representation of ,v, and is 
necessarily convergent. In fact, the proof of its convergence 
in every case must be included in the demonstration of the 
possibility of representing any function ,v by the series (3). 
Now the expression for v at time t, 1s ‘found by multiplying 
the respective terms of the second ‘members of (3), by the 


2 
quantities (a) ; 


a series which converges when ¢ is positive, (m,, m,, &c. 
being arranged in ascending order of magnitude). Hence, 
when ¢ is positive, the series for v is still more convergent 
than the series for ,v, and therefore (2) gives a convergent 
series for v for every positive value of ¢. If, however, the 
time considered be r before the initial instant, the series (4) 
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which is divergent. Hence it will depend on the degree of 
the convergence of the series for ,v, whether the expression 
of » in this case be convergent or divergent. In the latter 
case the distribution represented will be impossible, and 
therefore there will be no distribution from which the initial 
distribution can be derived, by spontaneous motion, in the 
time r. There are three cases of the initial distribution, 
which we must specially consider, 


1. If the convergence of the expression for ,v, or of the 
series of coefficients A A, 


be of a lower order than that of the series (a); that is, if the 
ultimate convergence of this series be less than that of (a), for 
every positive value of ¢, then for any finite value, however 
small, of +, the series 
Age’, 
will be ultimately divergent. Hence, in this case, for any 
time, however small, before the initial instant, the distribu- — 
tion will be impossible, and therefore the given distribution 
cannot be any stage but the first in a system of varying 

2. If the convergence of the series for ,v be ultimately of 


the sume order as that of the series (a); that is, if the coefh- 
cients A,, A,, &c. can be put under the forms 


where a,, a,, &c. form a series which ultimately has a con- 
vergence of a lower order than that of the series (a), according 
to the definition previously given; then, at atime 7’ before 
the initial instant, the distribution is represented by | 


v = La; cos (m,x + n,), 

which belongs to the first class considered, and is therefore 
an essentially primitive distribution. Hence, in this case, a 
finite age, r’ may be assigned to the given initial distribution. 
3. If the convergence of the series for ,v be of a higher 


order than that of the series (a); or if, for any finite value of 
t, however great, the series | 


2 2 


is ultimately convergent; then, for every finite negative value — 
of ¢, the series for v will converge, and therefore v will have 
a possible expression. Hence, in this case, no limit can be 
assigned to the age of the initial distribution. 
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_ The simplest case of distributions which belong to the third 
class, is that in which the expression for ,v is composed of a 
finite number of terms of the form eget , 


cos (mx +n); 


but we may also readily form series of an infinite number of 
terms, which will be sufficiently convergent to satisfy the 
condition stated for the third class, and therefore any initial 
distribution represented by such a series may be deduced 
from distributions existing previously for any length of time, 

These remarks are sufficient to indicate the nature of the 
question proposed. The details of the convergence or 
divergence of the series employed will depend on the values 
which must be assigned to m,, m,, &c. for satisfying the 
conditions of any particular problem. 

If, instead of the initial distribution of heat in the rod, from — 
—o to + being given, we have only a part, that for instance 
on the positive portion of the rod, we must, to make the pro- 
blem determinate, have some other condition given. In such 
cases, the part of the rod, over which the initial distribution 
is not-given may be conceived to be removed, and the second 
condition will then generally relate to the extremity of the 
part of the rod considered. In a previous paper (On the 
Linear Motion of Heat, vol. 111. p. 170), I have shewn how a 
solution of the problem may be obtained by determining the 
Initial distribution which must be made on the negative part 
of the rod, in order that the condition relative to the zero 
point may be fulfilled. Thus, if the temperature at the zero 
point, or the extremity of the part of the rod originally given, 
be constrained to be a given arbitrary function of the time, 
this constraint may be effected by producing the rod inde- 
_ finitely in the negative direction, and impressing on the part 
produced a certain initial distribution, determined by equa- 
tions (6) and (10) of the Article already referred to, in terms 
of the arbitrary function of the time, and the given initial — 
distribution on the positive side. In many cases, however, 
the distribution so determined will be impossible, either over 
the whole extent of the negative part of the rod, or from 
some point at a finite distance on the negative side, to an 
infinite distance. Notwithstanding this impossibility, the 
solution of the problem obtained by the method described 
above will still give a finite possible expression for the 
temperature of any point on the positive side, at any time 
subsequent to the initial instant, which will be the tempera- 
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ture that would be actually assumed by the point, if the 
temperature of the zero point had, by some external applica- 

tion, been constrained to satisfy the given condition ; and the 
- only:thing indicated by the impossibility of the distribution 
on the negative side will be that this constraint cannot be 
effected actually, by adding the negative part of the rod, with | 
a certain initial distribution of heat, to the given positive 
part. It is unnecessary to enter into this question separately 
here, as the details are very similar to those given above. It 
should be remarked however that the possibility or impos- 
pron | of the distribution on the negative part will depend 
entirely on the function of the time, which expresses the 
variable temperature of the zero point, and not at all on the 
given initial distribution on the positive part; and on this 
account, if the initial distribution on the negative part be 
impossible, all its subsequent forms will generally be impos- 
sible also. 

Before leaving this subject, it may be well to notice a 
point relative to the theory of isothermal surfaces, or surfaces 
of equilibrium, which involves considerations analogous to 
those with which we have been occupied above. | 
It is a known theorem that an attraction every where 
_ perpendicular to any given closed surface S, may be pro- 
duced by the distribution of a given quantity of matter m> 
over the surface, according to a law which is in every instance 
‘determinate ; but in general it will be impossible to produce 
the same effect by the distribution of matter over any surface 
in the interior of S, not coinciding with it. If, for instance, 
S' were the surface of a cube, or any surface containing points 
or edges, this would obviously be impossible, and it would — 
probably also be impossible in the case considered by Poisson, 
in which S$ is composed of two spherical surfaces.. In every 
‘ case there will be an infinite series of surfaces of equilibrium 
without S, becoming ultimately a series of spheres having the 
centre of gravity of m for their common centre, each of which 
is such that m exerts an attraction on any point in it in the 
direction of the normal. Hence, if one of these surfaces, s, 
be given, we may not only produce an attraction every 
where perpendicular to it, by matter distributed over itself, 
but by matter distributed over S, or over any surface be-— 
tween S and s, enclosing the former. Hence s is analogous 
to a distribution of heat which may be produced by a pre- 
viously existing distribution. But, unless S itself have the 
particular property we are considering in s, we cannot pro- 
duce an attraction perpendicular to s by any distribution on 
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a surface within S. Thus s is analogous to a distribution of 
heat which cannot be produced by any previous distribution 
existing at a time before it greater than a certain limit 7. 
Again, in some cases the series of surfaces of equilibrium, 
of which s is one, may be continued indefinitely inwards, till — 
we arrive at a surface enclosing no space, as for instance 
when s is an ellipsoid, in which case the surfaces of equili- 
brium are confocal ellipsoids, and the series within s may be . 
continued till we arrive at an elliptical disc. Such a case is 
analogous to that of a distribution of heat, which may be 
produced from a distribution existing an unlimited time 
before. . P. 


V.—ON THE INTENSITY OF LIGHT IN THE SHADOW OF A 
“VERY SMALL CIRCULAR DISK. 


Ir has been proved experimentally by M. Arago, that the 

central point of the shadow of a very small opaque circular 
disk is sensibly as bright as if the light were not intercepted. 
Fresnel has given a general explanation of this curious result _ 
in the Appendix to his Memoir on the Undulatory Theory of 
Light, (Mem. de l Institut., tom. v. p. 460: Ann. 1821-22), 
and states that Poisson has ‘proved the same thing by direct 
integration ; but he gives no hint as to where this proof can 

be found. 

The usual method of finding the illumination of a point 
appears in this instance to present some difficulty. © For if 
(fig. 1) A be the centre of the disk, Q the point in the shadow 
whose illumination is to be found, and we suppose the waves 
of light plane and parallel to the plane of the disk, also if 
PQ=4, it is easily seen =. the whole vibration at Q will 


be expressed by 
27 C [sin 


sin — (vt —~u).du. 


The superior limit being QD, a the inferior Qa, a being a 
point on the rim of the disk. 

The integration introduces at the superior limit the rather 
doubtful expression cos #: but, supposing this to vanish, as 
the nature of the case shews that it ought, for the vibration 
at Q cannot be affected by the small waves originating at very 
distant points of the general incident plane wave, still there 
is another difficulty ; for ze integration above indicated, if 


| erformed, ves 
P Ch (vt - Qa) 


for the total vibration at Q, = therefore the illumination 
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will be expressed by C*A’, which expression being inde- 

pendent of the radius of disk, is considered to prove the 
proposition. But the difficulty is, that this expression is true 

~ however large the disk may be, as well as when it is very 
small, which is manifestly absurd. It is true that the same © 
reasons which led us to neglect the quantity cos « would 


27 
teach us to neglect the expression cos + (vt — Qa), when the 


radius of the disk is very large, and so the illumination 
would become zero, as it ought: but this does not remove — 
the difficulty from the analysis ; for, according to the principle 
usually adopted, the illumination ought mot to vanish, when 
the radius of the disk is considerable ; because, since, when a 
general wave is broken up, the effect due to the small 
elementary wave is supposed to depend on the distance of 
the illuminated point from the origin of the small wave, and 
on that only, it follows that the effect of a small wave spread- 
ing from the rim of a disk whose radius is AP (fig. 2) in» 
exciting a vibration at the point Q, will be the same as that 
of a wave spreading from the rim of a disk whose radius 
is AP’ in exciting a vibration at Q'’, provided only that 
PQ= P'Q. In fact the common method does not take into 
account the effect of obliguity, and seems to involve the same 
error as if, in calculating the illumination of a point by a 
bright surface according to the principles of common Optics, | 
we were to omit the consideration of the fact that the quantity 
of light emanating in any direction varies as the sine of the 
angle of emanation. : 


This problem has been very elegantly treated in_ the 
Mathematical Journal (vol. 11. p. 141), and the result there 
obtained is the same as that found above, viz. that the illumi- 

_ nation is measured by C*A*: but the method there given 
seems to labour under precisely those difficulties which have 
been just alluded to. The fact seems to be, that the principle 
of Huyghens, which supposes the illumination of a point on 

- which a luminiferous undulation is incident, to be the same 
as that which would result from the combination of spherical 
wavelets diverging with equal intensity in all directions from 
the elements of the wave-front, cannot be looked upon as a 
physical principle, but only as an artifice rendered necessary 
by the state of analysis, and which will not always represent 
the physical conditions of the problem. 


In this paper I have adopted Huyghens’ principle so far 
_as this, that I consider the vibration of the ether at a given 
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point to be the sum of the vibrations which would exist if 
the wave-front were divided into elements and each of these 
elements were a centre of disturbance; but I do not suppose 
the wavelets to diverge spherically, but according to another 
law which [I shall investigate. 
With this view I shall first find an integral of the partial 
differential equation, which seems applicable to the case of 


light, in the particular case 7x which the vibration is sym- 


metrical with respect to a certain axis, and the feet rapidly 
diminishes as we recede from that axis. 


I shall assume, as the equation applicable to the vibrations 


of the ether, 
and let the axis of Z ‘be that about which every thing i is sym- 
metrical: let 2*+y°*=p, then ¢ is a function of z, p, and ¢; also — 


do an. dp 
dp) 


Now if @ diminishes very rapidly as p increases ; that is, if 
2 


9 
~ 


we may suppose p very small, the term p st will be neg-— 
ligeable as compared with or at all events we may 


p 
integrate upon that hypothesis, and if it should turn out that 
this is not the case, the method must be abandoned. 


2 
- Neglecting, then, the term p “s , our equation becomes 


= \ de + 4 7} (3). 
Let | o=usin =| (ot — 2), 


we have 


ty 

a 
| 

W 
| ] 
du\ 
— ut+4 —}, 
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and @ = ae a (et z).... (4). 


Let us now see whether p “f was rightly neglected: 


2 
2 


| 

will be small when p is small, proviged 
is not large: and this is a condition which is undoubtedly | 
fulfilled ; for, though A is very small, a is very large indeed, 
and v* also cannot be very different from a’. I think, there- 
fore, we are justified in considering that the expression (4) 
very nearly represents the true physical conditions of the 
problem: in order that it may do so, however, v must be 

greater than a, and for this I can give no physical reason ; 
I shall, however, assume that it is so, and proceed to apply 
formula (4) to the question in hand. 

If we take an annulus of the plane incident wave, whose 
radius is Vp, the area of the annulus will be dp, and if 


where 6 = the radius of the disk. And therefore the illumi- 
nation will be expressed by 


\! a 
a 2 2 
\v-a/) 


This expression diminishes very rapidly as 4 increases, 
but if 6b be small it approaches very nearly to the value 
2 2 

a” ( “ , which is its value if 5 = 0, or if the light was 
tT 
not intercepted. And this seems to agree sufficiently ac- 
curately with the experimental result. | 
_I will now apply the same method to express the illumina- 
tion at any point within the geometrical shadow. | 

If we take the projection of the given point on the disk 
as the origin, and call its distance from the centre g, and 


1 AQ (fig. 1) = f, the vibration at Q due to the waves spreading A 
= from the different points of this annulus, will be expressed by f 

mw? — a? : 
-— 

ral a? sin (ot —f) dp 4 

2 

b 
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measure 0 from the line joining this point and the centre, 
the equation to the circumference of the = will be 


sin’ + (r cos - gy = 
oF — 2rg cos = 
sin? 
p=g cos 20 + + 29 cos 6 — g’ sin’ 8). 


T° 2 2 


the illumination at the point in question I, it will be easily 


seen that 


where the limits of p are the above value and ow, and the 
limits of @ are 0 and z, if we double the integral, since only 
half the —_ 18 included between 0 and 7: : 


we-put », for shortness, and call | 


/ = en [92 cos 26+b2+2g cos 6 dé /(be-g2 sin? 6)] 
RJ 0 | 


The integration here indicated cannot be performed gene- 
rally, but there are some particular cases which are worth 
noticing. 


If g be so small that g and higher | powers may be neg- 
lected, the integral becomes 


= (1 2nbg cos 0) dO = =, 


or the illumination is the same as at the centre. 
If we take in the square of g, the integral becomes 


| cos 20 » vos @ dQ = (1 2nbg cos - COS 20 
0 *b’g’ cos’ 0) dé 

= (1+ J) | 

2.2 


| (14 nb q’). 


Hence, the centre is not so bright as an annulus taken very 
near to it. 


There is one other case which I shall consider, viz. for the | 
geometrical shadow, or when g = 0; in this case the integral - 


becomes 
(142 cos 20 -2nb® cos 20 
[ ag = om ag, 
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Shadow of very small Circular Disk. 
which is less than 7, because e” is <1, and 
[ cos a ~ 2nb* cos 20) d0=7, nearly, 
0 | 0 | 


since } is very small. Hence the line of the geometrical 
shadow will be less bright than the centre. Between these 
there must be at least one ring of maximum brightness, there 
may be more. | 
How far all this agrees with experiment I cannot say, but 
at least this method of considering the problem appears to 
embrace the two. great facts of the brightness of the centre, 


and the existence of rings in the shadow. ©. 


VI.—ON THE MOTION OF THE CENTRE OF GRAVITY OF 

BROKEN BODIES. | 

By W. Watton, M.A. Trinity College. 
THE fundamental principle adopted by Huyghens in the | 
determination of the Centre of Oscillation or Agitation: of 
Complex Pendulums, which forms the subject of the fourth 
part of the Horologium Oscillatorium, is enunciated in the 
following words: ‘Si pendulum e pluribus ponderibus com- 
positum, atque e quiete dimissum, partem quamcunque oscil- 
lationis integre confecerit, atque inde porro intelligantur 
pondera ejus singula, relicto communi vinculo, celeritates 
acquisitas sursum convertere, ac quousque possunt ascendere ; 
hoc facto, centrum gravitatis ex omnibus composite, ad 
eandem altitudinem reversum erit, quam ante inceptam os- 

cillationem obtinebat.” 

From the nature of the argument of this section of the 
Horologium Oscillatorium, it seems probable, although no~ 
where, I believe, expressly stated, and although in opposition 
to the direct meaning of the words which we have quoted, that 
Huyghens really supposed each of the disconnected molecules, - 
after arriving at its point of zero velocity, to be there perma- 
-hently retained; so that the ‘centrum gravitatis ex omnibus 
composite ’ should be considered as permanently assuming its 
original altitude from the date of the arrival of the molecule 
of greatest velocity at its maximum altitude. Such a modi- 
fication of the proposition is at any rate necessary to its truth. 
Che corrected form of the principle which we have mentioned 
Involves indirectly, and in a particular case, the principle of 
Vis Viva, a generalization which in fact resulted from the 
principle of Huyghens in the hands of John and Daniel 
Bernoulli. A demonstration of the principle of Huyghiens, 
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in its unmodified form, has been attempted by Professor 
_ Whewell, in a work entitled The Mechanics of Engineering, 
ip. 142, by the application of the general principle of Vis 
Viva. The reasoning however of the Professor is not — 
satisfactory, in consequence of his having omitted to take 
account of the idea of time, an essential element of the 
‘problem. 
The object of thie paper is to offer a few remarks on . the 
motion of the centre of gravity of the disconnected mole-. 
cules, supposing them not to be retained in their respective 
positions of zero velocity, but to. be capable of descending 
spontaneously from their positions of instantaneous rest. It 
is not necessary to suppose the body to have been oscillating, 
before the dissolution of the connections of its molecules, 
about a fixed axis in the manner of a pendulum, but we may 
imagine it to have descended from rest in any conceivable 
way; the remarks we are about to make bere generally 
applicable. | 
Let m, m',m".... represent the masses of the molecules, 
« “their initial distances below the horizontal line 
through the initial position (G) of their centre of gravity; 
xz, x,xz'....their distances below the same line at the mo- 
ment of the dissolution of their mutual connections ; then, if 
u,u,u.... denote the velocities of the molecules on the 
commencement of their independent motions, we have, by 
the principle of the Conservation of Vis Viva, 


92 


mu + mur 

= 2g {m (a —h) + m'(z' - +m" (a ~ 
but, by the property of the centre of gravity of any number 
of bodies, mh + 90, 


and 


where z denotes the of the centre of of 
the molecules at the moment of their mutual disconnection 


below the horizontal line through (G); hence 


12 112 


mu + + +....= 29 (m+m' +m" 


Again, let v, v,v'.... denote contemporaneous ‘velocities of 
the molecules in their independent movements, and y, y', ¥ 
......corresponding contemporaneous altitudes above the 
horizontal line through (#7); also let 4, k’, k’.... denote 
the values of y,y',y'.... at the instant of disconnection. 
We will commence with supposing the molecules to be 
constrained to move from the moment of detachment along 
smooth invariable lines, and not to be subject to mutual col- 


« 
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lision; then, by the theory of the motion of heavy particles, 

whether the paths of the molecules previous to disconnection 
be joined continuously with those which they afterwards 
pursue, or their motions be diverted at the instant of their — 
mutual detachment without loss of velocity imto different — 

directions, | | 

mv = mu’ — 2mg (y-k), = m'u? 2m'g (y' - 
__.we have therefore 


mv + m'v? + m'v? + + mu? + mut 
but, by the property of the centre of gravity of bodies, 
| mk + +m’'k' +....=0, | 
and my + m'y'+ +....=(m+m +m 
where ¥ denotes the altitude of the centre of gravity of the 
disconnected molecules at any time above the horizontal lin 
through hence | 
mo + | | 
= mu + mu? + mu? +... 29 (m+ m'+m' +...) (8). 
From this last equation it is evident that y can never have a 
value greater than that given by the equation 
from (1); 
and therefore it is impossible for y ever to exceed Z in mag- 
It remains to consider whether, or under what circum- 
stances, y will ever become equal to zx. Suppose y =2; 
then, from (3), we see that — | 
+ mv? + m'v’ +.... = 9, 
and therefore, contemporaneously, 
or, from the relations (2), : | | 
hence we see that, if ever y become equal to Z, every mole- 
cule must have at that instant a zero velocity, or be at the 
highest point of one extremity of its arc of oscillation. Sup- 
pose now that ¢, r, denote respectively, the time which m 


takes to move through the first and second portion of its arc of 
oscillation, the two portions being supposed to be separated 
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by the position of m at the instant of its detachment. Let 
(¢,7),(, +7 ) (t, + a represent like times in relation to the 
molecules m’, m’, m’' . Then, since v, v, v.... must be 
contemporancously equal to zero in order that y may be 
equal to z, the condition for such equality will evidently be 
expressed by a class of relations of the form 

tin(t+ +n (¢+7)=....(4), 
where n,n’, n’.... are all positive integers. The 
of satisfying these relations by any finite values of n, ',n'. 
as will be evident on a little reflection, will essentially depend | 
upon the forms of the paths described by the molecules sub- 
sequently to their disconnection. In case the molecules were 
to commence their motions at the moment of detachment 
along cycloidal arcs, starting at the lowest point of each, we 


should have 
and the relations (4) would become 


2n+1=2n'+1 
and therefore NM=N=N =... 


t=r=tH=re=t=r 


or the centre of gravity of the disconnected molecules would 
rise at each semi-oscillation of each separate molecule into 
a position on a level with (G). For the sake of another 
illustration, suppose that the molecules are perfectly elastic 
and keep rebounding above a horizontal : in vertical 
lines ; in this case we must put 7, 7, 7... .all equal to nero, 
and the relations will become. | 


(2+ 1)t=(7 +1) = 
or,sincee $$ u=gt, w=gt, u=gt.. 

(n+1)u=(n' + 1) = (n' + 1) 
if, then, w', uw’... . be quantities relatively commensurable, 
it is plain that these relations may be satisfied by admissible 
values of n,”',n’....; on the other hand, if the requisite 
commensurability do not subsist among the ‘projectile veloci- 
ties, it is equally evident that the relations cannot be satisfied 
by any finite values of n, n',n'.... 3 in the one case, then, © 
the centre of gravity of the disconnected molecules will from 
time to time, at finite intervals, arrive at its original altitude, 
while in the other case it will never attain to it in any finite 
time. It is needless to multiply examples of the impossibility 
of satisfying the formule (4) by any finite values of 2, n', n..- 
it being obviously rather a singular coincidence that the 
requisite commensurability should subsist without previous 
urrangement than a state of things of probable occurrence. 
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We have been supposing in the preceding observations 
that the disconnected molecules oscillate permanently in in- 


variable curves, and that they are free from a lability to 
‘mutual collision. We will now suppose the molecules, con- 
sidered as perfectly elastic, to be capable of assuming any 
irregular motions whatever, and of impinging against cach 


other. Under this more general consideration of the possible 
movements, the formula (2) will be no longer applicable, the 


values of v, v’,v.... obviously depending not only upon 


u,u,u.... and the altitudes of ascent, but likewise upon 
the effects of collision. ‘The formula (3) however, as we 
know by the application of the Principle of the Conservation 
of Vis Viva, still continues to hold good; we may there- 
fore shew, as before, that it will be impossible for the centre 
of gravity of the molecules to rise to its original altitude, 


unless the velocities v, v’, ¢ .... become simultaneously equal 


to zero. Let ¢, denote the interval from the instant of 


detachment to the arrival of the molecule m at its first 


position of zero velocity ; ¢,, ¢,,¢,....representing the inter- 


| 
vals between its successive arrivals at such positions: let 


t,,t/,t,....denote similar intervals in respect to the mole- 


cule m', and so on for the rest. ‘Then, in order that the 


centre of gravity may attain to its original altitude, we. must 
have | | 


the suffixes n, n', 2....denoting certain integral numbers. 


The probability that for any finite values of n,n',n'.... 
these indefinitely numerous relations should be verified, is 
evidently indefinitely small. We may therefore lay it down 


as ageneral practical principle, that if a mass of matter, either 


in its descent from rest or at the instant of its starting, be 


. dismembered into its constituent molecules, which subse- 


quently pursue paths prescribed by the accidental circum- | 


stances of gravity or collision, even supposing the molecules 
to be perfectly elastic and discarding all consideration of 
friction or atmospheric resistance, its centre of gravity will 
never again rise to a level with its original position. 

The truth of this practical principle, it will be evident, 


depends not upon any physical impossibility of the arrival of 


the centre of gravity at its original altitude, but upon the 
infinite mathematical improbability of the conspiracy of the 
necessary circumstances within any finite time. ‘The observa- 
tions which we have made above, in relation to the move- 
ments of the disconnected molecules of descending bodies, will 
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evidently be equally applicable to the case of fluids flowing 


down from a higher into a lower level, the centre of gravity 
of which will accordingly oscillate in a lower basin without 


_ ever reaching its original altitude. 


If instead of supposing the mass to be dismembered into its 
constituent molecules, we imagine it to be broken merely into 
fragments of finite magnitude, since the equation (3) will still 
hold good, we shall still have relations such as those of (5) 
for the conditions of the rise of the centre of gravity to its 
original height ; it is important to remark, however, that the 


number of these relations will, according to the new hypo- 


thesis, be finite ; for if three molecules of each fragment, 
which are not in a single straight line, ever assume zero 


velocities, the whole of each fragment will be in a position of 


instantaneous rest. ‘Thus, if there be p fragments, there will 
be only 3n—1 relations of condition; if the mass be supposed 


~ never to have fallen to pieces at all, then, according to the 


same reasoning, we see that there will be only two relations 
of condition. ‘Thus we come to this general conclusion, that 
if a mass of matter in descending from any position fall to 
pieces, the improbability of its centre of gravity ever rising 


to its original altitude will be the greater, the greater be the 


number of the fragments ; and that, when the fragments 
degenerate into molecules, the improbability will amount to 


- practical impossibility. It may perhaps not be superfluous 


to observe that, if all the particles of a descending body be 
constrained to move in parallel planes, the whole body will 
be at rest if any two molecules, not in the same horizontal 
line, be at rest, and therefore the equations of condition will 
be reduced to one. If the body oscillate about a fixed axis 


as a complex pendulum, then if any one molecule be at rest 


the whole body will be for a moment without motion; and 
therefore there will be no relations of condition, or the centre 
of gravity will undoubtedly rise at the end of each oscillation 
to its original altitude. 


VII.—ON THE INVERSE CALCULUS OF DEFINITE INTEGRALS. 
By GeorGE 


THis paper must be considered as the continuation of a 


paper on the Transformation of Definite Integrals, (Journal, 
No. XVII. p. 216.) 


M. Liouville has given the following theorem as the basis _ 
of an inverse calculus of definite integrals, viz. 
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a demonstration of which will be found in an able article in 
this Journal, vol. 1. p. 118. It is there observed, that if the 
form of ¢(z) be such as to render the first side infinite, the 
second may be made to agree with it by the aid of the com- 
plementary function, so that the theorem is true for all forms 
of p(x). M. Liouville, I believe, considers his result as only 
true under a limitation, and this appears to me to be the 
correct view. ‘There are at least two distinct cases of the 
theorem, each of which admits of a convenient transformation | 
in the case of m being fractional. 
Let us consider the integral {dz z""g(z +a), assuming in 
the first instance that @(x) may be developed in ascending 
powers of z, and that the first exponent exceeds -1. 
Put xz = - az, and let a = then 


[de (a + 2) =(-fa" -2)} 
= (1 -z)} 
= {dz z)% (6°). 


If we take for the limits of z, 0 and 1, the corresponding 
limits of z will be 0 and - a, and the above gives 


J 


—+n+4 1 

a result unquestionably true for all positive values of ». Let 
_n be an integer, then | 


= (-)" eal (a) 


(a+ 2) = (27). .. .(29), 


= (-)"P(n) . (30), 


a conclusion easily verified by integration, it being observed 
that the limits of the integrations relative to @ are 0 and a. 
We might obviously regard this theorem as true for fractional 
values of m; but, without thus extending our assumptions, we 
may see from the above, that M. Liouville’s formula does 
| G2 
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not apply when ¢(z) is only susceptible of development in 
ascending positive powers of | 

Let us now suppose the development of $(z) a descending 
one, and the exponents negative. 


Put z =az, as before, 

[dx (a+ x) [dz 2"" (1 + ...(81), 
wherein 5 is of negative interpretation. Now, by a well- — 
known theorem, (Gregory’s Examples, p. 473), | 

| dzz*" (1 + 2)" = (32), 


provided that n be positive and less than 7. Hence, com- 
paring with the second member of (31) taken between the 
same limits, 


whence, observing that the limits of z are also 0 and x, | 


\ 


Hence, it appears that M. Liouville’s formula applies when 
(x) involves only negative exponents in its development, 
whereof the first with its sign*changed exceeds n. ‘The 
limits of the integrations relative to a in the second member 
are+manda. § 6. 

Finally, let us suppose x frational. 

The object of investigation being to determine the function 
under the sign of integration, when the definite integral 
is given, let us assume that within the supposed limits, 
[dz x" = a known function, and that we wish 
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8 to ascertain the form of (a). In the case of the ascending 
development we have from ( 29), on changing a into ©, 


| 
| d 
dé 


Let 7 be the integer next below x, being 0 when n is a 
proper fraction, 


r 1): 
(2 ) n+1) 
qi"! 


T(t-n+1) 
=(-) | ave (1-v)'"Y(av)...(37) 


Similarly for the descending development of ¢(z), and for 
negative exponents, we have, by (33), 


 (<€) 
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P(e) = 
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Now let z be the integer next below , then 


d0\ d0 


| -—4+n-7)Te-n+1) 
: d d 1 ( = 
=(-)" —— ant — 


I'(n) da da 
| 


We shall somewhat simplify the above formula by changing’ 
2 into 2— 1, and in (30) and (37), x into — z, in order to avoid 


the negative limit. Our results will then be comprised in the 
following theorems. | ; 
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1. If, being a positive fraction, | “de g(a - x)= (a), 


0 
and if 7 being the integer next above n, 2 - 1 be less than the 


first exponent in the ascending development of ~(a), then — 


i-n 


but if 2 be a positive integer less than the same first exponent, 


1 d n | 
(a) = P(a). 


2-4, If, » being a positive fraction, [ dzxz"'o(a+xj=W(a), 


the exponents. in the development of ‘L( a) being all negative, 


and if 2 be the integer next above m, then 


but if » be a positive integer, | | 
The application of these formule, as it presents no difh- 
culties, so it requires no illustration. . They enable us, when- 
ever (a) satisfies the necessary conditions, to determine (qa) 
either in the form of a definite integral, or in that of a series ; 
but for the latter object it is perhaps simpler, changing z — 
into <’, to employ at once the untransformed equations (35) 
and (38). 
I apprehend that the same principles may be extended, by 
aid of the theorem of M. Dirichlet, to the inverse calculus of 
definite multiple integrals. At present, however, I have not 


leisure to pursue this inquiry, nor am I so much concerned 


to multiply results as to establish principles. 
Lincoln, Oct. 26, 1842. 


~ 


VIII.—ON A NEW SPECIES OF EQUATIONS OF. DIFFERENCES. 


THE sort of equation to which the heading refers is that in 
which the degree of the equation with respect to one variable 
is defined by the other, as in wu, =u, oe where wv, is a 
function of the integers z and y. A problem of some interest 
leads to such an equation, and will tend to shew the sort of 
solution which may be given. ) 

It is well known that the number of ways in which the 
number » may be put together is 2”": thus 4 can be con- 
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structed in 2° ways; namely, 4, 31, 13, 22, 211, 121, 112, 
1111. But in this theorem orders are counted; thus, 112 is 
considered as different from 121, and both from 211. ‘The 
corresponding problem, namely, to find in how many ways 
the number z can be constructed, different orders not being 
considered as different ways, is of much greater difficulty, 
and [ am not aware. of any thing having been written 
upon it. 
In the first sins, 3 if & be not less than n, the number of 
ways in which z can be put. together is the same as that in 
which » + 4 can be composed of numbers one of which is /: 
thus 16 can be 8+...,. in as many ways as 8 can be put 
together. If then uw, represent the number of ways in 
which z can be made up of y, and numbers not exceeding y, 
the problem is solved when w,, , is found: for the number of 
_ ways in which z can be put together is u,, od whatever | 
~ may be the value of 4, from 0 inclusive. 
Again, it is obvious that the number of ways in which an 
number, as 16,;can be composed of numbers no higher than 
- 6 and not all lower, or u,,.,, 18 the sum of the ways of com- 
posing 10 out of numbers not exceeding 1, 2...6 Every way 
in which 10 can be so formed, gives 16 by adding 6. Hence 


Z.y x 


y. 1 z-y.y? 
Hence 
OF 


-an equation of the y- degree with respect to 2. 


This equation gives a tolerably easy mode of constructing 
a table: for all numbers up to 10 it is as follows; wu, ,, wu, ,, 
and , being always units. 
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The rule being, after the units are written down, from the 
y'" number of the 2" line by adding the (y - 1) number of 
the (z- 1)" line to the (x — y)* number in the z line. Thus 
10 can be formed out of numbers one of which is 4, and no 
other exceeds 4, in 9 ways: these 9 ways are 442, 4411, 433, 
4321, 43111, 4222, 42211, 421111, 4111111. - 

As soon as we can express u,,_,., In terms of z and y, we 
have then an equation of the y order for determining w, 
Say that wu =V_,, we have then, 


y-l 


y 


which, considering y as a constant, is of the form 
w.—w. =V_. 
| x 


From this consideration we can arrive at the general form o 
the solution, and this is all I shall here attempt. First we 
have = 1, so that 


2 U2, 2 1; 


the complete integral of which, the constants being found 


from vu, ,= 0, uw, ,=1, is 


the complete integral of which, determining the constants 


from 0, .= 0, 1, 
6z°- 7-9(- 1) + y’) 
79 ’ 


where (3 and y are the imaginary cube roots of unity. This 
. May be expressed as follows:—The number of ways in 
_ which z can be composed of numbers, none greater and not 
all less than 3, is | | 

C £31 Cet — 1 

12 12 12 12 12 12 
according as x divided by 6 leaves a remainder 0, 1, 2, 3, 4, 
or 5. If we were now to attempt to generalise, we might 
suppose that | | 


u,., is of the form A,, + P,+....+P,, 


where A represents a rational and integral function of the 
n” degree, and P. is a periodic or circulating function of 
the z'" order, or going through a cycle of » values, But this 
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inference would not be quite correct. The terms of each 
successive integration are of the form 3 


definite constant x (u*V_), (u" = 1). 


_ When PV, is a periodic term of the form Cj3*, (3 being some 
previous root of unity (say the m'), generally speaking 

is a simple circulating function of the m™ kind ; 
aad if Vv. be of the form A 3’, the form resulting from 
integration is also 43". But i B = uw, which happens among 
the cases arising when n is a multiple of m, then pd." V, 
takes the form 4" A, or the function of z is raised a unit in 


dimension. Hence the most general form of the solution of | 


=A, + + + 


whine a, the of the of ; Is the 
integer in n - y divided by y. 


The next step gives for w,. , the denominator 864, and the 
following numerator : | 


62° + 182° — 27x — 39 + (27x + 27) (- 1y 


this is the 144" part of 


(6) 2° + 3x 36, 
+ 82° - 9274 5,- (1) + 82? = + 5, 
(2) + —-20, (8) 2° + 32° + 96, 

(3) 2° + 82° - 9x - 27, (9) a + 32° - 9x - 27, 
(4) 2° + 82° + 32, (10) + 32° - 4, 


(5) + 92-11, or (11) 2° + 82" — Oz — 11. 
according as x divided by 12 gives the remainder 0, 1, 2.... 
or 11. 

It thus appears: that u.,, may be described as the in- 
Leger nearest to z + 12, and u.., as the integer nearest to 
(2° + 34°) + 144 or (2° + - 9x) 144, according as x 1S 


even or odd. Probably this simple species: of description 
“might be continued. 


-IX.—NOTES ON MAGNETISM. NO. lI. 
By R. L. Eviis, M.A. Fellow of Trinity College. 


A GEOMETRICAL construction, by means of which the action 
of a small magnet on a distant particle of free magnetism may 
be readily determined, is mentioned in a memoir by Weber, | 
on the Bifilar Magnetometer (Screntific Memorrs, 11. p. 270). 
~Itis due to Gauss, but I do not know where he hae demon- 
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strated it. A proof of it may be acceptable to some readers 
of the Journal. . 

I begin by enunciating the construction in question, which 
will be easily understood without a figure. Let AB be a 
small bar magnet, ¢ its centre, P the particle of magnetism 
on which it acts. Join cP, draw vid perpendicular to it, 
meeting eP produced in D. Let cQ=}cD. Join PQ. Then 
_ PQ or QP (according to the sign of the magnetism of P and 
the direction of the poles of AB) is the direction of the action 


of AB on P; and = - is its ‘magnitude, M being the 
measure of the magnetism of 4B, m that of the magnetism 
of P. 

The dimensions of the magnet being small, and its length . 
in the direction of its axis being much greater than its 
breadth or thickness, we may proceed as follows. 

Conceive the magnet to be composed of a series of intense 
particles ranged along its axis. Let s be the distance of any 
one of them from ec, uds the measure of its magnetism. Also 
let cP. =r, and call the angle cP makes with cD, 0. 

Then. the distance of the particle in question from r,s 

(r° + s° — 2rs cos 8), and consequently its action on P is 


muds 
r+ 2rs cos 
magnetic attraction being supposed to follow the ordinary 
law. The component of this action along cD, is 
muds | 


(7 cos - s). 
8? - 2rs cos 
This is 3 approximately equal to 


or to met fr cos + (3 cos’ 1) s}ds; 
r 
and therefore the total action of AB on P, parallel to cD, 1s 
| 8 cos’ 1 


the integrals being taken wn the whole so a of the 
magnet. Consequently = 


since the aggregate magnetism of a magnet, each element 
being taken with its proper sign, is zero. Again, /usds 
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- being the moment of the magnetism of AB, is the measure of 
its magnetic power, or what we have called M. 
Consequently the action parallel to cD is 


The action of the element at s, perpendicular to cD, 1s 
muds | 


- r-sin 
(7° + — 2rs cos 0) 
or, approximately, 


(1 + 3 - cos 8) sin 6ds. 
r r 


The total action perpendicular to cD is, therefore, 


sin cos 8: ..... 


The equation to the resultant of these two forces is, since the 
line — through P, | 
y —rsin 
3 sin 0cos@ 
For y = 0, or at the point Q, 


or sec 0. 
Now cD =r sec 9, therefore 
cQ=1cD, 
which proves the first part of the construction. 


Next, to find the magnitude of the whole action, square 
and add (1) and (2), then — 


(1 - 6 cos’ 8+ 9 cos* 8+ 9 sin’ cos’ 
| 
Mm 
0) 


is the magnitude sought. Now 
PQ =r {sin’ 0+ (cos } sec 9)*}4 
r(1-—243 sec’ 0) 
1 (sec? + 3), 
and cQ=!rsec 
+ 3.cos 0). 


‘Therefore 
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Consequently, if R be the magnitude of the resultant sought, 
Mm PR 


tr” OF 
oo which was to be proved. 


The preceding formule enable us to determine all the 
circumstances of the mutual action of two magnets, which are 
such as to fulfil the conditions of our hypothesis. 

- For instance, in the memoir Intensitas vis Magnetice 
Terrestris, Gauss has shown that if a magnet and a needle 
be placed at right angles to one another, then the moment of 


rotation of the needle due to the action of the magnet is 


approximately twice as great when the line of the axis of the 
magnet passes through the centre of the needle, as it is when 


the line of the axis of the needle passes through the centre of » 


the magnet. 

‘In neither case has the transverse force (2) any tendency to 
produce rotation. Its effect is destroyed by the fixed centre 
of the needle. 

Let y be the distance of any element of the needle from the 
centre, mdy its magnetism. ‘Then, in the first case, we shall 


have @ = A approximately, (R being the distance between 


the centre of the magnet and that of the needle). Conse- 
quently the force (1) may be expressed by the formula — 


Mmdy (2 - y \ 


since cos’ = 1 — nearly; 


Rs) 
which, neglecting y’, becomes 
2Mmdy 
the moment of this round the ceutre of the needle is a dy a 


and the total moment is, therefore, 
2MM' 
In the second case, = 0 for every element of the needle, and 
the moment sought is, therefore, | 
MM. Mm 
| 3 


since 
b 


and the sign is immaterial. ‘The moment in this case is, there- 
fore, one half of what it was before, which was to be proved. 


where M’ = {mydy. 


following general investiga-_ 


tion, in which we shall ascertain the snient of rotation rie 
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to the action of one magnet or needle upon another, whatever 
be their relative positions. 

The data by which we shall suppose the position of the 
magnets to be determined, are the distance between their 
centres, the angles which the axis of each respectively makes 
with the line joining their centres, and the angle between the 

axes themselves. The last element may be readily replaced 
_ by the dihedral angle between two planes which intersect in 
the line joining the centres of the magnets, and one of which 
passes through the axis of one of the magnets while the 
other passes through that of the other. — 

Let the axis of the magnet, whose action on the other is to 
be calculated, be taken as axis of z, the centre of the magnet 
being the origin of co- -ordinates. Let x, y, z be the co- ordi- 
nates of any point in the other magnet. ‘Then, if r*=27+y'+2’, 


+ 2° 
3 cost @-1 = and sin cos @ = 


Again, let a, b,c be the co-ordinates of the centre of the 
second magnet, a, (3, y the cos angles its axis makes with the 
co-ordinate axes, p the distance of any element dp from the 
centre, X, Y, Z the forces on dp parallel to the axes of 
co-ordinates, m the intensity of dp. ‘Then 


= (82° - dp, 
Y= xydp, 
Z= xzdp 


Let G, H, K be the moments of these forces about linie 


drawn through the centre of the second magnet parallel to 
the axes of co-ordinates, 


H= (x - a) - (82° (z dp, 
(y b) - szy(z- a)} d 
Now t=a+ap, y=b+p, z=C+y¥p. 
Hence, neglecting the ‘square, &c. of p, 
G= a (by cf3) 
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Notes on Magnetism. 
H = {3aca A y} pdp, 
{(8a* — B - 3aba} pdp, 


Integrating for p, we find, for the total moments, 
83MM MM' 


(K) = === a (af - ba) 


MM. being for the second magnet what MM is for the first. 
Let L be the eg of these three moments, then 


+ (H) + (KY. 


Consequently, a’ + (3° + 1, we shall have 


(aa (aa + UB + oy) 


Now, let 0, 6 be the angles which the axes of the two mag- 
nets make respectively with the line joining their centres, 
and let @ be the angle between the axes themselves. Then 


a= Roos @ and aa+bB+cy = R cos also a =Ccos @. 


Consequently 
{9 cos*6 sin’'+ 6 cos (cos cos + cos 0) + sin*9 |, 
or L= + 3 cos —(3 cos 8 cos — cos 


Let x be the dihedral angle already mentioned, then 
cos @ = cos 8 cos 0 + sin @ sin 9 cos x. 
Consequently the last equation becomes 


{1 +3 cos’ (2 cos 6 cos - sin sin cos x) }4, 


the required expression. 


? 
Bs 
me 

if 

an 

Sor 

a 

an 

We 

/ 
ws 
is 

oad 

ae x, 

4 

a 

5 

ie 

‘ ‘ 

/ 2 

\ 
4 \ 

\ 

ay 

a 

he 

= 
a 
L 
£4 
. 


ADDENDUM TO ART. Il. 


The method of this paper may be applied to the integration 
_ of any linear equation with constant coefhcients, by means of 
the following theorems, the first of which answers to integra- 
tion by parts, the second to John Bernoulli’s Theorem. Let 
D denote the operation of finding a differential coefficient, 

then 
(D+ ay. PQ= P(D + a)'Q- (D+ {P'(D+a)'Q}, 
P being the result of DP: and, abbreviating D + a into 9, 
0". PQ = Pe'Q- Pe’*Q+ 

which, if P be rational and integral, gives 0” PQ: in terms 
of O 1, &e. 

Anne, to find all the preceding equations of — 

(D-1)(D- 2fy=X, | 

those of the 13" degree, for instance: operate upon 

(D- 1)" - 1)°(D - 2)'y} =(D - 1)" X), 

(D 2)" {a" (D 2)°(D - = (D- 2)" (2"X), 
for all integer of maid n, from m= 0 to m= 5, and 
from x= 0ton=7. We shall thus have the 14 equations of 
the 13 degree ; ua from each of these, exactly as in the 

preceding part of this paper, others of the 12", &c. degrees 
may be obtained. 

The criterion which, being satisfied, shows that 

= {P (D+ 

will when operated upon by the symbol (D + a)", not show 
that symbol in the result, precisely resembles:the correspond- 
_ Ing criterion when « = 0, or the case of simple integration. 

The preceding method may be equally applied to finding 
the intermediate equations of linear partial differential equa- 
tions with constant coefficients, or of equations of differences. 

With respect to equations of differences, the process 
answering to integration by parts is as follows. Let £ 


signify the operation of changing z into «+ Az, then, P 
and Q peg functions of z, 


(E + a)'.(PQ) = AP.(E+ a)' EQ}, 
the extension of John Bernoulli’s theorem is, E+ 4 
being 9, | | 
0'(PQ) = AP. 6° EQ + A’Pe*E’Q - 
which terminates when / is rational and integral. 
ih applying this to such an equation as 


- ay" - by'....y + X, 
instead of using 2* as a ellintier, it would be more con- 
venient to use the factorial get”, A. DE M. 
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